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Abstract —This paper considers proportional fairness amongst 
ACs in an EDCA WLAN for provision of distinct QoS require¬ 
ments and priority parameters. A detailed theoretical analysis is 
provided to derive the optimal station attempt prohahility which 
leads to a proportional fair allocation of station throughputs. The 
desirable fairness can he achieved using a centralised adaptive 
control approach. This approach is based on multivariable state- 
space control theory and uses the Linear Quadratic Integral 
(LQI) controller to periodically update CWmin till the optimal 
fair point of operation. Performance evaluation demonstrates 
that the control approach has high accuracy performance and 
fast convergence speed for general network scenarios. To our 
knowledge this might be the first time that a closed-loop control 
system is designed for EDCA WLANs to achieve proportional 
fairness. 


1. Introduction 

Enhanced Distributed Channel Access (EDCA) was pro¬ 
posed in IEEE 802.1 le-2005 standard to support QoS enhance¬ 
ment and service differentiation for WLAN applications [1]. It 
extends the basic Distributed Coordination Eunction (DCE) by 
classifying traffic flows into four different Access Categories 
(AC), namely voice, video, best-effort and background. Traffic 
with higher QoS requirements, e.g. shorter delay deadline, is 
assigned a higher priority, and hence, on average, waits for 
less time before being sent to the channel. This mechanism 
is beneficial for high-priority traffic. Compared to the DCE, 
EDCA sacrifices the performance of low-priority traffic to 
some extent to provide QoS support for high-priority traffic. 
When the network is saturated with a large proportion of high- 
priority flows, an extremely unfair scenario will appear, in 
which the channel will be almost completely occupied by high- 
priority flows, e.g. VoIP or video streaming flows, however 
low-priority traffic, such as email or web browsing data, will 
suffer severe starvation. 

Resource allocation in EDCA WLANs has therefore been 
the subject of considerable interest. The objective is to seek 
for a fair allocation of network resources (e.g. throughput, 
airtime and etc.) amongst different traffic types, and mean¬ 
while, guarantee the specific QoS requirements and service 
differentiation. This paper considers proportional fair allo¬ 
cation of station throughputs amongst ACs for provision of 
distinct average delay deadlines and priority parameters. The 
802.1 le EDCA standard specifies four contention parameters 
to distinguish priority levels, which are minimum Contention 
Window {CWmin), maximum Contention Window {CWmax), 
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Arbitration Inter Erame Space (AIFS) and maximum Trans¬ 
mission Opportunity (TXOP). A set of default values for the 
four parameters are recommended in the standard for each 
physical (PHY) layer supported by 802.lie. As the default 
values do not take into account the varying WLAN conditions, 
and thus lead to suboptimal performance and no fairness 
guarantees, in this paper we find the optimal CWmin value 
that leads to proportional fair allocation of station throughputs 
while assuming AIFS and TXOP taking the recommended 
values and CWmax = CWmin- The optimal CWmin value 
corresponds to an optimal station attempt probability which is 
derived from the proportional fairness analysis. 

In order to implement the derived proportional fair alloca¬ 
tion in practice, a centralised adaptive approach which uses 
multivariable state-space control theory is then proposed. The 
WLAN is represented as a discrete multi-input multi-output 
(MIMO) linear time-invariant (LTI) state-space model. A state 
feedback control method, the Linear Quadratic Integral (LQI) 
control, is used to tune the CWmin value to drive the station 
attempt probability to the optimum so as to maintain a fair 
throughput allocation. We have demonstrated in simulations 
that the proposed control approach is adaptive to general 
network scenarios with high accuracy and fast convergence 
speed. To our knowledge this might be the first time that a 
closed-loop control system is designed for EDCA WLANs to 
achieve proportional fairness amongst ACs. 

The remainder of this paper is organised as follows. Sec¬ 
tion II gives a comprehensive review of the state-of-the-art 
research on fairness and control theory approaches to solve 
network problems. Section III presents the theoretical analysis 
to derive the optimal station attempt probability which leads 
to proportional fair allocation of station throughputs given the 
constraints on the average delay deadlines for different ACs. 
Section IV describes a centralised adaptive control approach 
which can realise the proportional fair allocation derived from 
Section III in real networks. Section V evaluates the perfor¬ 
mances of the fairness algorithm and the proposed centralised 
control approach, and Section VI concludes the paper. 

II. State of the art 

A. Fairness 

Eairness has been the subject of a considerable body of lit¬ 
erature on 802.11 WLANs [2]-[19]. The unfairness behaviors 
may be caused by a number of factors, e.g. hidden termi¬ 
nals, exposed terminals, capture, uplink/downlink unfairness, 
asymmetric radio conditions and multiple data rates and etc.. 
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which have been investigated in [2], [4], [5], [6], [8], [7], [11], 

[13] . There also exist distinct fairness criteria that are widely 
adopted in network resource allocation, such as time-based 
fairness, throughput-based fairness, proportional fairness, max- 
min fairness, weighted fairness and etc. (see for example [11], 
[18], [17], [19], [10], [9], [13], [14] and references therein). 
In this paper, we employ the proportional fairness criterion 
to deal with the unfairness amongst ACs with differentiated 
priority parameters. 

The CSMA/CA scheduling used in 802.11 differs funda¬ 
mentally from wired networks due to carrier sense deferral 
of the contention window countdown and the occurrence of 
colliding transmissions, both of which act to couple together 
the scheduling of station transmissions and lead to the rate re¬ 
gion being nonconvex [20]. Therefore, well established utility 
fairness techniques from wired and TDMA networks cannot 
be directly applied to random access CSMA/CA wireless 
networks, and hence most studies on proportional fairness in 
802.11 WLANs are confined to approximation approaches [2], 
[3], [13]. [18] corrects prior approximate studies and provides 
the first rigorous analysis of proportional fairness in 802.11 
WLANs. It shows that there exists a unique proportional fair 
rate allocation and completely characterises the allocation in 
terms of the total air-time quantity. [17] extends the work in 
[18] by considering lossy links and BSC-based coding with 
delay deadline constraints. The optimal joint allocation of air¬ 
time and coding rate allows the throughput/loss/delay tradeoff 
amongst flows sharing network resources to be performed in 
a principled manner. The proportional fairness analysis in this 
paper builds upon the approaches used in [17] and [21]. 

In particular, the fairness issue in 802.lie EDCA WLANs 
has been given considerable attention. [11] derives a through¬ 
put allocation based on the proportional fair criterion in 
multirate 802.1 le WLANs. It shows that in a proportional fair 
allocation high and low bit rate stations are assigned with the 
same share of channel time, and thus high bit rate stations 
obtain higher throughput. Two schemes that respectively in¬ 
volve Contention Window CW and Transmission Length TL 
(which is based on TXOP) are then proposed to achieve 
this allocation. [13] investigates the weighted proportional 
fairness in both single-rate and multi-rate 802.1 le WLANs 
via test-bed experiments, and compares proportional fairness 
with time-based fairness in a multi-rate setting. Proportional 
fairness can be achieved by tailoring the CWmin, and the 
time-based fairness can be achieved by adjusting packet size or 
TXOP limit. It concludes that in a multi-rate 802.1 le WLAN 
proportional fairness with equal weights achieves higher per¬ 
formance than time-based fairness in terms of both aggregate 
utility and throughput. [11] and [13] deal with the unfairness 
behaviour arising in a WLAN due to asymmetric channel 
conditions, and [13] investigates the problem through test-bed 
evaluations relying on no theoretical analysis foundation. [12], 

[14] , [16] address the unfairness problem existing amongst 
ACs. The priority-based service supported by 802.lie EDCA, 
while allowing differentiated service for flows of different 
priorities, cannot ensure service amount in proportion to their 
demands. [12] proposes a mechanism called Weighted Eair- 
EDCA (WE-EDCA) which uses Distributed Eair Scheduling 


(DES) in each backoff entity to provide weighted proportional 
fair service among different ACs. [14] proposes an algorithm 
to compute the optimal configuration of the EDCA differenti¬ 
ation parameters given a set of QoS requirements in terms of 
throughput and delay with multiple real-time and data ACs. A 
throughput and delay analysis is provided, based upon which 
the optimal configuration algorithm is derived by maximising 
the throughput using the weighted max-min fairness criterion. 
[16] presents a new scheme which exploits differentiations 
of both inter frame space and contention window to achieve 
weighted fairness for two classes of services under EDCA 
mode in an 802.1 le WLAN. Given the AIESs, the proposed 
scheme can properly set the corresponding CWs such that 
the ratio of the two classes’ successful transmission proba¬ 
bilities can attain a pre-defined weighted-fairness goal. [15] 
proposes a dynamic contention window control scheme to 
achieve fairness between uplink and downlink TCP flows for 
the IEEE 802.1 le EDCA-based WLANs while guaranteeing 
QoS requirements for real-time traffic. The proposed scheme 
first determines the minimum contention window size in the 
best-effort access category at APs, and then determines the 
minimum and maximum contention window sizes in higher 
priority access categories, such as voice and video, so as 
to guarantee QoS requirements for these real-time traffic. In 
this paper, we provide a comprehensive throughput and delay 
analysis that incorporates EDCA differentiation parameters for 
802.1 le EDCA WLANs and derive a throughput allocation to 
achieve proportional fairness amongst ACs. 

B. Control theory approach 

Control theory has been applied to the area of commu¬ 
nication networks in a wide range of aspects. Eor instance, 
[?]nalyses a combined TCP and Active Queue Management 
(AQM) model from a control theoretic standpoint. It uses 
a nonlinear dynamic model of TCP to design a feedback 
control system depiction of AQM using the random early 
detection (RED) scheme. [23] introduces a control theoretical 
analysis of the closed-loop congestion control problem in 
packet networks. The control theoretical approach is used in 
a proportional rate controller, where packets are admitted into 
the network in accordance with network buffer occupancy. A 
Smith Predictor is used to deal with large propagation delays. 
[24] proposes a QoS-provisioning feedback control framework 
in order to achieve TCP uplink/downlink fairness and service 
differentiation. The medium access price (MAP) is delivered 
to TCP senders and the TCP senders adjust their sending 
rates to reduce congestion at the interface queue of the home 
gateway in an 802.11-based home network. [25] proposes a 
centralised adaptive control (CAC) approach to dynamically 
adjust the CWmin configuration of 802.11 WLANs with the 
goal of minimising the overall throughput performance. A 
proportional integrator (PI) controller is used to establish a 
closed-loop control system. [26] extends the work in [25] by 
considering real-time traffic in 802.1 le WLANs. The CWmin 
configuration is adjusted in terms of minimising the average 
delay, which results in a better quality of experience (QoE) of 
the video traffic. [27] proposes a distributed adaptive control 
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Fig. 1. Network model. 


(DAC) algorithm based on the multivariable control theory. 
An independent PI controller is installed at each station in 
a WLAN and uses locally available information to drive the 
overall network performance to the optimum. In addition to 
the closed-loop control methods, [29] proposes an open-loop 
self-adaptive rate control approach for multi-priority WLANs. 
The approach is based on a biological competitive model 
which guides data flows to compete for network bandwidth 
in the way of a native ecosystem. The model parameters 
self-tune themselves to optimise the bandwidth utilization 
in an EDCA WLAN with multiple ACs. In this paper we 
employ the multivariable feedback control method to design a 
centralised adaptive control approach to achieve proportional 
fairness amongst different ACs. The WLAN is represented as 
a discrete MIMO LTI state-space model. The state feedback 
LQI controller is used to tune the CWmin value. To the best 
of our knowledge, this might be the first time that a closed- 
loop control system is used in 802.11 WLANs to deal with 
fairness issue. 


TABLE I 
Notation 


N 

Number of ACs 

rii 

Number of STAs in AC i 

n 

Total number of STAs 

Wi 

Minimum contention window of AC i 


Maximum contention window of AC i 

AIFSi 

Duration of AIFS of AC i 

SIFS 

Duration of SIFS 

EIFS 

Duration of EIFS 

Trts 

Duration of RTS 

Tcts 

Duration of CTS 

a 

Duration of a physical time slot 

u 

Number of time slots in AIFSi 

tmin 

Minimum t value among all ACs 

mi 

Number of packets in the TXOP burst of AC i 

L 

Packet size 

r 

PHY data rate 

n 

Station attempt probability of AC i 

pidle 

Slot idle probability 

psucc 

i 

Probability of a successful transmission by a station of AC i 

psucc 

Probability of a successful transmission in a time slot 

Si 

Station throughput of AC i 

pcol 

Duration of a collision 

psucc 

Duration of a successful transmission from a station of AC i 

i 

Protocol overhead of a TXOP burst of AC i 

poo 

Protocol overhead of a single packet within TXOP burst 

M 

Retry limit in exponential backoff algorithm 

pcol 

Conditional collision probability of AC i 

pblk 
^ i 

Blocking probability of AC i 

Df 

Expected countdown delay of AC i 

phlk 

Expected blocking delay of AC i 

pretx 

Expected retransmission delay of AC i 

ppucc 

i 

Expected retransmission delay of AC i 

Di 

Average delay of a TXOP burst of AC i 

di 

Delay deadline of a single packet of AC i 

H 

Plant lumbarisation matrix 

x{k) 

System state at instant k 

y(fc) 

System output at instant k 

u(k) 

System input at instant k 

r(k) 

Controller input at instant k 

K 

LQI optimal gain matrix 

Q 

LQI state cost weighting matrix 

R 

LQI control cost weighting matrix 


III. Proportional fairness in multi-priority 
802.1 lE WLANs 

A. Network model 

We consider a single-hop 802.1 le EDCA WLAN with one 
AP and n client stations, as depicted in Pig. 1. The channel is 
assumed to be error-free for all supported PHY rates. Traffic 
flows are classified into N different ACs. We assume that each 
client carries flows of a single AC, so there are no virtual 
collisions in our setup. The number of stations in the ith AC 
is Hi. The total number of stations is thus n = The 

analysis can be readily generalised to encompass situations 
where client stations have lossy links and carry more than one 
AC. But the current simplified error-free model is sufficient 
to capture performance features of differentiation settings in 
WLANs. 

B. Station throughput 

We start with the analysis of station throughput under 
saturation conditions, i.e. there is always a packet waiting to 
be transmitted in the queue of a station. Por the ith AC, the 
following parameters are defined; is the minimum 


contention window; is the maximum contention win¬ 

dow; ti is the number of time slots in the period of AIFSi, 
i.e. AIFSi = SIFS + ti x a where SIFS is the duration of 
the Short Interframe Space and a represents the duration of a 
physical time slot; rrii is the number of packets transmitted in a 
TXOP burst. We assume that packets of all ACs have the same 
length of L bits, and are transmitted at the same PHY rate r 
Mbps under the assumption of error-free channels. Due to the 
use of TXOP bursting, the RTS/CTS exchange mechanism is 
used to make fast recovery from collisions. The notation used 
in this paper is listed in Table. 1. 

A MAC time slot may either be a PHY idle slot, a successful 
transmission or a colliding transmission. Let Ti denote the 
probability that a station carrying a flow of AC i attempts 
to transmit in a time slot, so 0 < Ti < 1. The probability that 
a time slot is idle is 

N-l 

P^dle = [| (1 - Ti)"* 

i=0 

As the channel is assumed to be error-free, packet losses are 
only caused by collisions. The probability that a station with 
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a flow of AC i makes a successful transmission is then 


P- 


= n{i - - 


N-l 


■ n 

3=0,j^i 


A jjidle 

1 - A 


The probability that a time slot is a successful transmission is 
then 


N-l 


N-l 


psucc _ 


riiTi 


2=0 


psucc _ pidle \ ^ ‘ 

'H^i — ^ / > -I 


2=0 


The throughput of a station carrying a flow of AC i is thus 
given by 

psucc^ . T 

s-(t) - _^ ^ _ 

- N-l 

pidleg- _j_ psucc'psucc -|- (^1 pidle psucc'^'pcol 

2 = 0 


in which = T^xs + PIPS is the duration of a collision. 
EIFS represents the duration of the Extended Interframe 
Space used in 802.11 WLANs, which is given by EIFS = 
Tack + SIFS + DIFS. Note that T™' is the same for all 
ACs due to the use of RTS/CTS handshaking. is the 

duration of a successful transmission from a station sending 
traffic of AC i. It depends on the size of TXOP packet burst, 
and is thus given by 

J.SUCC = T° + 

where T° = Trts + SIFS + Tqts + AIFSi is the protocol 
overhead associated with the transmission of a TXOP burst; 
T°° is the protocol overhead associated with each packet 
transmission within a TXOP burst, i.e. T°° = TpHYhdr + 
2SIFS + Tack. 

By working in terms of the quantity ai = , cti > 0 

instead of r^, the station throughput is rewritten as 

_ aimiL 

^ . J'col ' ^ 


in which 


N-l 


N-l 


iV J . rriQiinc 


i=0 


i=0 


C. Average delay 

Next, we will calculate the average delay experienced by 
a TXOP burst of each AC. We start with the analysis for 
ordinary 802.11 MAC scheduling with binary exponential 
backoff algorithm and then move onto the scenario when 

= CW^ax. 

The delay is defined in this work as the duration since a 
station starts contending for the medium until the transmission 
is finished (either received successfully or dropped because of 
reaching the maximum retry limit). The calculation is based 
on the EDCA WLAN throughput model derived in [28]. The 
average delay consists of four expected delays, described as 
follows: 


• Expected countdown delay. Eor each backoff stage g 
S 9 S M, and M is the retry limit. We assume that 
CW^ax ^ the average countdown delay 

for AC i is in which CWi,g = 


is the contention window at the gth backoff stage. The 
expected delay associated with the backoff countdown 
process is then given by 

/ ^ 9 r'lAi 

Df = - pr') E 

^ g—O h—Q 

/2 = 0 ^ 

where is the conditional collision probability for the 
ith AC in the throughput model, i.e. 

iV-l 

pcoi = 1 _ (1 _ p[ (1 _ (2) 

3=0 

Expected blocking delay. During the countdown process, 
when a transmission is detected on the channel the 
backoff time counter is “frozen”, and reactivated again 
after the channel is sensed idle for a certain period. A 
station is called blocked in our analysis when it senses 
(an) ongoing transmission(s) from some other station(s) 
during its countdown process. The blocking delay is the 
period during which a station is “frozen”. The expected 
number of time slots in the backoff countdown process 
is jcr. At each time slot, a station could be blocked 
by either a successful transmission or a collision. Eor 
a station of AC z, the delay caused by a successful 
transmission from some other station is 


N-l 






N-l 


N-l 


E Tr^ngTgil - J] (1 - - a)”*-' 


3=0 


fc=0 

k¥=3,i 


The blocking delay because of a collision is 

. N-l 

Dbc ^ J.COI h _ (1 _ 7-^)"^-! (1 - r,)"^ 

^ 3=0 

3¥=i 

Af-1 

- (n* - l)Ti(l - n 

3=0 

N-l N-l 

- E n - Ti) 

j—O k=0 


The expected blocking delay of AC z is thus 

jjcd 

J^blk ^ _^i^jgbs 

a 

Expected retransmission delay. The expected retrans¬ 
mission delay for AC z is calculated by multiplying 
the expected number of retransmission attempts by the 
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collision duration, i.e. 

j-^retx _ rpcol ^ 

. M s 

Y,g{pr'y{i - pr‘) + {m+ i){prT^^ 

\g =0 ^ 

• Expected successful transmission delay. The expected 
successful transmission delay is the duration of a suc¬ 
cessful transmission multiplied by the probability that the 
transmission is not dropped, which is given by 

JJSUCC _ rJ^SUCC^ _ 

Combining the above four delays, the average delay of a 
TXOP burst of AC i is therefore given by 


The proposed approach in this paper works by finding the 
optimal contention window to achieve proportional fairness 
amongst ACs , so the exponential backoff algorithm is unnec¬ 
essary in our setting and we simply set 
i.e. M = 0. To simplify notations, we hereafter refer to 
Wi. The four expected delays then become: 


A-= 3 




d: 


_ rpCOL pt 


col rycol 


DSUCC ^ _ p. 


<col J 


According to the throughput model in [28], when 
CWmax = CWmin the Station attempt probability under 
saturation conditions can be reduced to 

^ 2{1-Pn . 3 . 

* 2(1 - Pf'=) + W^-1 ^ ^ 

in which 


pblk _ 




N-l 


ti 




]=o 


is the probability that the backoff counter is suspended due 
to a busy channel during the period of AIFSi. tmin is the 
minimum t value among all ACs. 

Similarly by working in terms of the quantity ai = jzfy, 
the average delay experienced by a TXOP burst of AC i when 
M = 0 is then 


+ 


_ ic*(cr-f 
“ 2 

WrY, / 

2(l + a,)”‘-^'^ ’ 


_ ^ _ /pSUCC pCol\ _|_ pcol 

' ( 1 -f * 

- -f - r“')(n, - l)a*) 

(4) 


in which 


N-l 


Y^= + 


N-l 

Z,= Y. {Tr^-T^°')n,a, 

and 

Wi =—(p + ai) '^{l + aj) - 1-1 

1=0 

D. Proportional fair allocation 

The 802.lie EDCA standard provides service differentiation 
by assigning different contention parameters to distinct ACs. 
Delay-sensitive traffic flows, such as voice over WLANs and 
streaming multimedia, are assigned with higher priorities. 
This mechanism has a significant cost for lower priority 
traffic flows as they can practically starve in dense network 
deployment. In this section we aim at finding the optimal 
a := [ai]jg{o,i, " ,Af-i} to achieve fair allocation of station 
throughputs amongst ACs. Meanwhile we take into account 
the delay constraints for each AC. The utility function is 
defined as the sum of the log of station throughputs 

Af-l 

U{a.) := ^ nilogSi(o) 

i^O 

Di{oL) < rriidi 0 < i < N — 1, 
ai > 0 0 < i < N — 1. 

in which the station throughput is given by Eqn. (1); the 
average delay is given by Eqn. (4); di is the delay deadline 
for a single packet in a TXOP burst of AC i. 

By plugging in the station throughput expression and remov¬ 
ing the constant terms, the optimisation problem is simplified 
as 

N-l 

max Ufa) := 7 ni(log cti — log X) 

i^O 

s. t. Dfa) < niidi 0 < i < N — 1, 
ai > 0 0 < i < N — 1. 

1) Non-convexity: It can be verified by inspection of the 
second derivative that the objective function is not concave 
in a and hence the maximisation problem is not a standard 
convex optimisation task. We proceed by making the log 
transformation rji = loga^. The optimisation problem then 
becomes 

max Ui{rf) := Ufe^) 

^ (5) 

s. t. Dfe^) < rriidi 0 < i < X — 1 . 

Lemma 1. Ufri) is concave in r). 

The proof of Lemma 1 is included in the Appendix. 

2 ) Solving the optimisation with KKT conditions: To solve 
this problem, we will use KarushKuhnTucker (KKT) condi¬ 
tions. The Lagrangian is 

N-l 

L = C/i(77) - ^ (A(e^) - mid^) 

2=0 


max 

OL 

s. t. 
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where the Lagrange multiplier fXi > 0. Differentiating the 
Lagrangian with respect to rji and setting it equal to zero yields 




n dX dDi 

X dr]i diji 


,N-i} 

in which 


Af-l 

j=0 


dD^ 

% 


- n* = 0 


(6) 


dr]i 


'Y'SUCC 

^ i 

J^col 


X N-1 

- Ijrtio;* +ni (1 + aj)"" 

^ 3=0 


Oti 

1 -|- OCi 


% V 2(l + ai)"*-i 


+ 


((T + T“')aA dW^ , 


dai 


(1 + ai) ”• {ui - 


1 ) 


(^(jsucc _ - n^)a, + 1 ) - Z, + 

+ Y,a,{Tr^ - T“')(l - n,)(l + a,)"”* 


dDj _ - r“0(^i - l)aj + % - T™') 


i9?7i 


2(1 + 


((T + T“')a*A 51L, 


9a, 


+ a,(l + aj)^-”^- (r/““ - T“‘) 


.IL, 


IL, 


col 


dY^ 

9 a, 


ILllLai si_ndZj 

The derivatives 4^, and are respectively given 

aaj ’ aoij ’ aa^ acnj •' ® 

by 


9 VL, 

9ai 



N-l 


1 + ai) (1 + ctj) 


3=0 


rii — 1 . , \ 

Z ; (f, imin “t“ 1 ) 1 

1 + ai / 


ti 


-i 




+ 1 


dW, 

9 a i 


2?7/' / 

= ■ a,(l + a,)^^^' ■ 

^-1 f -f 4-1 

n \ I'j l^min I J- 

(1 + afe)-"'') 


91^^ 

9ai 

9ai 


ti(l + ai) 


^rpSUCC _ rj-iCol'^ 


3) Subgradient algorithm for optimal a: Given the values 
of the Lagrange multipliers p*, the solution to Eqn. ( 6 ) speci¬ 
fies the optimal a. To complete the solution to the optimisation 
it therefore remains to calculate the optimal multipliers p*. 
These cannot be obtained in closed form since their values 
reflect the network topology. We proceed in a centralised 
manner by using a standard sub-gradient approach. The dual 
problem for the primal problem defined in Eqn. (5) is given 
by 

min g{p) 

fi>0 


where the dual function g{fi) is given by 


7V-1 


g{p) = max U'{ot) + V' gi{m,di 

CX. 

2=0 

A ^-1 

= f/^(a*(/x)) + ^ fiiirriidi 


2 = 0 


Eor any a, 


Di{a)) 

A(a*(M))) 


N-l 

g{p) > U'{oi) + y] lii{midi - Di{oL)) 

2 = 0 

and in particular, the dual function is larger than that for a. = 
i.e. 


N-l 

9{^) > U'{oL'{fjL)) + ^ gii{m^di - Di{oL'{fji))) 

2=0 

A ^-1 

= g{fJi) + ^ (/X 2 - - Di{oL'{fi))) 

2 = 0 

A sub-gradient of g{-) at any fi is thus given by the vector 


midi 


A(a*(A)) 


ie{0,l,-- ,N-1} 


and the projected sub-gradient descent update is 




- 7 • (niidi - Di{a* 


where 7 > 0 is a sufficiently small stepsize, and [/(-j]^ := 
max{/(-),0} ensures that the Lagrange multiplier never goes 
negative [31]. 

The subgradient updates for p can be carried out centrally 
by the AP. Eor the zth AC, the AP requires the knowledge of 
the number of stations rii, the PHY date rate r and the packet 
size L. The algorithm to calculate optimal a is detailed in 
Algorithm 1. 


Algorithm 1 Calculate optimal a 

1 : Initialise t = 0. 

2 : repeat 

3: f = f-h 1 ; 7(9 = 1/^2 

4: solve for a*(/x(9j by combining equations 


/,(a,/r(9)=0 VzG{ 0 ,l,--- ,iV-l} 

5: Vt G {0,1, • • • , W — 1}, calculate 




(t+i) 


1 " - 7'‘> 


{rriidi - Di[a* 


6 : 


until 


(i) G ,iV- l},\mid^ - < e, 

where e > 0 and is sufficiently small; 

(ii) Vf G {0,1,-- - ,N - l},mid^ - D^(a* > 0 . 


IV. Centralised closed-loop control approach 

In this section, we design a centralized adaptive control 
approach to implement the desirable proportional fairness in 
real networks. Based upon the analysis in Section III, the 
proportional fairness is achieved when the station attempt 
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probability parameter a. reaches its optimum value a*. The 
variable a is only determined by the minimum contention 
window Wi with AIFS and TXOP taking the recommended 
values and CWmax = CWmin- Our approach uses a multivari¬ 
able closed-loop control system to tune W to drive the station 
attempt probability to its optimum. As the station attempt 
probability is hard to measure in real networks, we measure 
the conditional collision probability instead of ol{t) 

in the proposed control approach. 


A. Measuring p^°’' 

The measuring of p‘^°^ is performed periodically at the 
AP every beacon interval. It requires messages passing from 
ordinary stations. A transmission from a station collides if 
a CTS is not received within a predetermined period after 
sending a RTS. Each station starts to count the number of 
transmitted RTS packets and that of missing CTS packets 
after receiving a beacon packet. The numbers counted are 
piggybacked to the AP along with following RTS packets. 
After a beacon interval, based on the piggybacked information, 
the AP knows the number of transmitted RTS packets from 
stations of AC z, denoted as and the number of 

missing CTS packets of AC i, denoted as , within 

a beacon interval. If we ignore the transmission duration of 
the beacon packet, and assume that packets collide with a 
constant collision probability within each interval, the average 
number of collided packets of AC i is These 

packets will be observed as missing CTSs, which yields 

. Therefore, we have 


E 


i^M-CTS\ 
N^RTS j 


jyRTS Ri 


The observed conditional collision probability p°^^ can be 

pj-M-CTS 

estimated as 'pfUTs ■ To simplify notations, we hereafter 
refer to p'^°^ with p. 


B. Control algorithm 

The closed-loop control system consists of two modules as 
depicted in Fig. 2. 

The controller module is installed at the AP and periodically 
carries out the adaptive algorithm every beacon interval, which 
is typically 100ms in 802.11 WLANs. The adaptive control 
algorithm is described in Algorithm 2; 

The plant is the WLAN itself. The input of the plant 
is the contention window W = [Wqi • ’ ’ and the 

output is the observed conditional collision probability = 
■ ■ ^Pn-i]- The design objective is to obtain a stable 
system in closed-loop with desired performances and shape 
the output of the system to the given reference value. The 
reference value is the optimal conditional collision probability 
p* given by Eqn. (2) for t — t*. 

C. Linearisation of the non-linear plant 

As the proposed adaptive control algorithm is executed 
every beacon interval, the period is long enough to assume that 


Algorithm 2 Adaptive control algorithm 

1 : AP calculates the optimal a* and the resulting p*. 

2: AP broadcasts the current W to all subscribed stations 
along with a beacon packet. 

3: Each station starts to count the number of transmitted 
RTS packets and the number of missing CTS packets after 
receiving a beacon. The numbers are piggybacked to the 
AP along with following RTS packets. 

4: At the end of a beacon interval (100ms), AP counts the 
total number of transmitted RTS packets of AC i, 
and the total number of missing CTS packets, _ 

The observed conditional collision probability p°^^ in this 
interval is then calculated based on the counted numbers. 

5: The reference p* and the observed are input into the 
controller to calculate a new set of W. 

6 : If the output of the controller W is not an integer, it is 
rounded to the closest integer value, and it has to be at 
least larger than 1. This requires hardware modihcation to 
allow CW to be chosen as an arbitrary integer. 

7: Go back to step 2. 


P‘ 


Controller 


W 


WLAN 


pobs 


Fig. 2. Closed-loop control system 


the measurement corresponds to stationary conditions. This 
implies that depends only on the current W, i.e. the 

system has no memory. Following this. 


N-l 

pfs = !_(!_ (1 _ ^ 7 ) 

7=0 

3=^1- 

in which is a function of Wi, given by Eqn. (3). 

Eqn. (7) and Eqn. (3) give a non-linear relationship between 
pobs pj/- jjj order to simplify the controller design, we 
proceed by working with a linear approximation to this non¬ 
linear relationship around the stable point of operation. 

The perturbations of input around the stable point of oper¬ 
ation is 


w = w* + sw 


in which W* is the W which yields the optimal value r* 
from Eqn. (3). 

The perturbations suffered by p°^® can be approximated by 
Sp°^^ =SW H 
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in which 



dpt‘ 

^p'n-i 

awo 

awo 

dWo 


apf‘ 

^p'n-i 

awi 

awi 

awi 

9pn'’= 

dpf‘ 

9PnLi 

OWm-i 

awN-i 

aiVjv-i 


The partial derivatives can be respectively calculated as 

^ dpf^ dTk 

dW, ~ ^ dTk ' 
k—0 

and 

^ ^ 9pf" dTk 

^ ~ to ^ ^ 


in which 


dpObs 

dn 


N-1 


n (1 - 

fc =0 


rii -1 


dpObs 

dTj 


N-1 


11(1 


Tii 


k^O 


{1-Ti){l-Tj) 


The control task can be accomplished by using the LQI 
control method [30] to design our controller. Fig. 3 shows the 
control block diagram for the system, in which x{k) G is 
the system state, y{k) € is the system output, u{k) € 
is the controller output and r{k) G is the controller input, 
which is the optimal collision probability p*(k). K € 
is the control gain matrix, and B G and C G 

are the state-space system matrices. Tg is the sampling period 
of the system, i.e. the beacon interval 100ms. 

The LQI controller computes an optimal state-feedback 
control law by minimising the quadratic cost function 

OO 

j{u{k)) = Y, {k)Qz{k) + vF {k)Ru{k)) 
k=0 

for any initial state a:(0), in which z{k) = [x{k);s{k)] and 
s{k) = s{k — 1) -f Ts • (r(fc — 1) — y{k — 1)) is the output of 
a discrete integrator. 

The matrices Q and R are the weighting matrices respec¬ 
tively indicating the state and control cost penalties. Q and R 
are required to be real symmetric and positive definite. 

The state feedback control law is defined as 


dW, -2(1 - (1 + (n, - 1)(L - + l)n) 

dn _ _ Tj{l - Tj) _ 

dWj -2(1 - Pf^)ni{l - Tj){tj - tmin + 1) 

At the stable point of operation, = r*, the non-linear plant 
is thus linearised as 

p°bs ■H{t*)+p* ( 8 ) 


u{k) = —Kz{k) 

The optimal state feedback gain matrix K is computed by 
solving the associated discrete algebraic Riccati equation 

P = A^{P - PB{R + B'^PB)-^B'^P)A + Q 

in which 

A = [0jVx2Af; —C * Tg JaTxTv] 


D. State feedback control 

With the lumbarisation, the WLAN can be represented as a 
discrete MIMO LTI state-space model. According to the pro¬ 
posed adaptive algorithm, the conditional collision probability 
at instant fc -f 1 is determined by the contention window input 
to the WLAN at instant k, the state and measurement equations 
are therefore given by 

J x{k -f 1) = Bu{k) 

\ y{k) = Cx{k) 

in which the system state is the conditional collision proba¬ 
bility, 

x{k) = [p°^%k)]^ 

the system input is the minimum contention window, 

u{k) = [Wik)]^ 

and the system model matrices are 
B = -H^ 


and 


C — Inxn 


The system output is thus 

y{k) = Cx{k) = [p°^%k)f 


and 

B = [B; O^TxAf] 

K is constructed from the solution of the above algebraic 
Riccati equation P* and weighting matrices Q and R, which 
is given by 


K = {R + B^ P* B)-^ B^ P* A 


E. Selection of Q and R 

The selection of weighting matrices Q and R affect the 
performance of the LQI controller. A simplified form using 
only 3 degrees of freedom is chosen in this work. The matrices 
are of the form 


and 


Q 


QiIn 

0 


0 


92 In 


2Nx2N 


(9) 


R = p ■ Inxn 


( 10 ) 


in which qi G is the weight for the state feedback cost; 
<72 e R+ is the weight for the integral feedback cost and 
pGR+ is the weight for the input cost. 
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Fig. 3. LQI Controller 


TABLE II 

802.1 1 PROTOCOL PARAMETERS USED IN THE SIMULATIONS 


cr 

9 /is 

TpHYhdr 

20 /is 

SIFS 

16 /is 

TrtS 

46.67/ls 

DIFS 

34 /is 

Tcts 

38.67 /LS 

EIFS 

88.67 /LS 

Tack 

38.67 /LS 


TABLE III 

EDCA CHANNEL CONTENTION PARAMETERS FOR 802.1 1 OFDM PHY 


Access categories 

AIFSN 

Max TXOP 

Background (AC_BK) 

7 

0 

Best Effort (AC BE) 

3 

0 

Video (AC_VI) 

2 

3.008ms 

Voice (AC_VO) 

2 

1.504ms 


V. Performance evaluation 
A. Throughput and delay performance 

The main objective of this work is to achieve proportional 
fair allocation of station throughputs while satisfying specific 
delay constraints of different ACs. To verify if the proposed 
fairness algorithm meets this objective, we first evaluate the 
throughput allocation and delay performance. The results are 
obtained using Matlab based on the throughput and delay 
analysis in Section III-B, III-C and the proportional fairness 
algorithm described in Section III-D. The throughput and delay 
analysis is based on the 802.1 le performance model presented 
in [28] under the assumptions that stations have saturated 
traffic and CWmax = CWmin- The accuracy of this network 
model has been fully verified in [28] for different network 
scenarios. 

An example of an 802.1 le WLAN with traffic of two ACs 
is considered, one of which is data traffic belonging to AC_BE 
(best effort), and the other is video traffic belonging to AC_VI 
(video). Flows of both ACs are saturated. The 802.11 OFDM 
PHY layer is assumed to be used. The PHY data rates for 
two ACs are the same, i.e. r = 54Mbps. The packet size 
is L = 8000 bits. The 802.11 protocol parameters used in 
the evaluation are listed in Table II. The FDCA contention 
parameters recommended for 802.11 OFDM PHY layer are 
listed in Table III. Note that the IFFF 802.1 le standard 
also provides recommended values for CWmin and CWmax- 
As our control approach searches for the optimal contention 
window to achieve proportional fairness, the default CWmin 
and CWmax are not used in the proposed approach. 

Fig. 4(a) and Fig. 4(b) shows the throughput and delay 
performance for two ACs versus the number of stations in 
AC_VI while keeping the number of stations in AC_BF fixed 
as 1. We let the TXOP burst reach the maximum limit as listed 


in Table III. The average delay deadline for a single video 
packet is d 2 = 250/is, which is the successful transmission 
duration of a video packet. The delay deadline for a single data 
packet is four times that of a video packet, i.e. di = 1000/rs. 
It can be seen that under these delay deadline constraints, the 
resource allocation can be divided into four phases. Phase I; 
When there are one video station and one data station, the 
network load is quite light, and both delay deadlines are not 
reached yet. Phase II: When the number of video stations 
increases to 2, the increased collision possibility leads to 
longer delays, and so the delay deadlines of both ACs are 
reached. As the number of video stations increases, in order 
to achieve a fair throughput allocation, data stations attempt 
to access to the channel more with an increased attempt 
probability, while video stations attempt less with a slightly 
decreasing attempt probability. The fairness algorithm makes 
the throughput of two ACs get closer to each other. Phase III: 
When the number of video stations increases up to 5, it comes 
to the turning point when video traffic is so aggressive that the 
proportional fairness algorithm allocates higher throughput to 
video traffic by assigning increased attempt probabilities to 
both ACs, but the increase for video traffic is larger than that 
for data traffic. The delay of video traffic is then reduced to 
be less than the deadline limit, while the delay constraint of 
data traffic remains tight. Phase IV: As the number of video 
stations continues increasing, the throughput ratio between 
video and data traffic remains around 1.5 in this phase. Fven 
with 10 video stations and only one data station, the data 
station can still deliver a reasonable amount of throughput. 
Fig. 4(c) shows the corresponding station attempt probabilities. 
It can be seen that contrary to the 802.1 le FDCA standard, our 
algorithm assigns data traffic with a higher attempt probability 
although it has lower priority. However not only does the delay 
performance satisfy the QoS requirement, the throughput is 
also fairly allocated between two ACs. Fig. 4(d) shows the 
corresponding collision probabilities for both traffic types. As 
expected, the increased number of video stations results in an 
increase in collision probability in general. In Phase IV the 
collision probability for data traffic increases slowly, while 
that for video traffic slightly decreases. 

B. Air-time 

The presence of collision losses and the coupling of station 
transmissions via carrier sense make the flow air-time in a 
WLAN not simply be the successful transmission duration but 
also include air-time expended in collisions. We define the^ow 
total air-time as the fraction of time used for transmitting a 
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(a) Throughput 



(b) Average delay 



(c) Station attempt probability 



(d) Collision probability 

Fig. 4. Throughput and delay performances in an 802.lie WLAN with two 
ACs, the number of stations in AC_BE n\ = 1, r = 54Mbps, L = SOOObits, 
di = 1000/is and d2 = 250//S. 


TABLE IV 

Comparison of flow total air-time allocations under 

DIFFERENT DELAY DEADLINE CONSTRAINTS. 


Flow 

1 

2 and 3 

4 and 5 

6 

Sum 



VI VO 

BK 


Case I 

0.1565 

0.1530 

0.1550 

0.1562 

0.9287 

(!^ase II 

oTSST 

0.1667' 

0.1667 

0l667^ 

1 


flow, including both successful transmissions and collisions. 
For a flow of AC i, the flow total air-time is 

TfSUCCrpsUCC i ncolrpcol 

rpair _ _ ^2 ^i _ ~r ^ _ 

^ “ N-l 

pidle^ _|_ psuccpsucc -|- _ pidle _ psucc'jpcol 

i^O 

1 / / , Ti \ 

pidle ) 

The work in [17], [18] finds that the proportional fair 
allocation assigns equal total air-time to each flow in a 
WLAN, and the air-times sum to unity. We investigate the 
air-time allocation by considering an 802.1 le WLAN with 
four ACs, in which AC_VI and AC_VO have two stations, 
i.e. n 2 = ns = 2, and AC_BE and AC_BK have one station, 
i.e. ni = 714 = 1- The PHY rates for four ACs are the 
same, i.e. r — 54Mbps. The packet size is L = 8000 bits. 
Table IV compares the flow total air-time allocations under 
different delay deadline constraints. In Case I, the average 
per-packet delay deadline for the four ACs are respectively 
di = 900/7S, d2 = 300/rs, ds = 250/rs and di = 1800/rs, 
while in Case II, the average delay deadlines are relaxed as 
di = ds = ds = d^ = 5000/is. 

It can be seen that different from observations in [17], [18], 
the flow total air-time is not equalised in Case I, and the sum 
of air-times is less than 1. We also notice that in this case the 
delay constraints for voice traffic Flow 2 and Flow 3 are tight. 
Nevertheless when the delay constraints are relaxed in Case II, 
none of the delay deadline constraints is tight. Flows are now 
allocated with equal total air-times, and the air-times sum to 
1 . It is thus found that the air-time allocation in our algorithm 
is affected by the imposed delay constraints. With tight delay 
constraints, the proportional fair allocation assigns flows with 
the exact amount of air-time needed by each of them. Air¬ 
time resource in the network is not completely occupied in 
such a case. However with loose delay constraints, flows can 
occupy all the available air-time resource, and air-time is 
evenly distributed amongst flows in a network as discovered 
in previous work. It is worth pointing out that since the flow 
air-time usage overlaps due to collisions, the flow total air¬ 
times summing to unity does not imply that the channel idle 
probability = 0. 

C. Q and R tuning 

The tuning of qi, <72 and p can be performed using the trial 
and error method. The influences of the three parameters are 
illustrated using an example with three ACs, AC_BE, AC_VI 
and AC_VO. The number of stations in AC_BE, AC_VI and 
AC_VO is respectively ni = l,7i2 = 2 and ns = 1. Three 
ACs use the same PHY rate, i.e. ri = r 2 = rs = 54Mbps. The 
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packet size is I = 8000 bits. The average packet delay limit 
for AC_BE, AC_VI and AC_VO is respectively di = 900^s, 
^2 = 300/rs and = 250/rs. The TXOP burst reaches the 
maximum limit. Fig. 5 plots the system output response for 
AC_BE with different sets of qi, q 2 and p values. The reason 
we do not put the output responses of AC_VI and AC_VO 
in Fig. 5 is that we notice that three outputs have the same 
convergence speed with a fixed set of qi, q 2 and p values. The 
effects of the three parameters are outlined as follows: 

• qi imposes the constraints to the state dynamics. It is 
directly related to the overshoot. A higher qi corresponds 
to a lower overshoot. As shown in Fig: 5(a), qi = 700 
results in an overshoot, while qi = 800 corresponds to 
an undershoot. 

• q 2 impacts on integral action dynamics and so on the 
system dynamics. As shown in Fig: 5(b), the higher it 
is, the smaller the rising time will be, and the higher the 
overshoot will be. 

• p affects the dynamics of the controller input, and so 
on the system dynamics. It is related to the overshoot. 
A higher p results in a higher overshoot, as shown in 
Fig: 5(c). 

D. Adaptivity to changes in the WLAN 

We will next evaluate the adaptivity of the proposed method 
to changes in the network size. The scenario being considered 
is depicted in Fig. 6. The algorithm starts at f = Os with four 
saturated stations, one in AC_BE, two in AC_VI and one in 
AC_VO. One more station in AC_VO joins the network at 
t = 100s and leaves at f = 200s. At t = 300s one AC_BK 
station joins the network, and after 100s one AC_VI station 
leaves. The PHY data rates for the three ACs are the same, i.e. 
r = 54Mbps. The packet size is I = 8000 bits. The average 
packet delay limit for data, video, voice and background traffic 
are respectively di = 900/is, d 2 — 300. /rs, ds = 250/is 
and di — 1800/rs. Fig. 7 plots the variation of contention 
windows over time. Fig. 8 plots the corresponding station 
throughput for each AC. Q and R take the form as displayed 
in Eqn. (9) and (10). We choose qi = 750, q 2 = 2000 and 
p = 0.005 to make a fast convergence speed. It can be seen that 
when the network condition changes the contention window 
converges to the desirable value very quickly as long as proper 
Q and R are chosen. Moreover, the steady-state errors can be 
neglected, which means the control system has high accuracy 
performance. 



(a) 92 = 2000, p = 0.005 



(b) 91 = 750, p = 0.005 



(c) 91 = 750,92 = 2000 

Fig. 5. System output step response, n\ = 1, n 2 = 2, 723 = 1, r = 54Mbps, 
L = 8000 bits, di = OOO/j-s, d 2 = SOO/rs dz = 250/2S. 


VI. Conclusions 


This paper considers using a closed-loop control approach 
to achieve proportional fair allocation of station throughputs 
in a multi-priority EDCA WLAN. The optimal station attempt 
probability that leads to proportional fairness is derived given 
the average delay deadline constraints of different ACs present 
in an WLAN. To achieve the desirable proportional fairness, 
a centralised adaptive control approach is proposed. The 
WLAN is represented as a discrete MIMO LTI state-space 
model. The LQI control is used to tune the CWmin value to 
the optimum. We have demonstrated using numerical results 


1 AC_BE Sta 1 Ac_VO Sta 1 AC_VO Sla 1 aC_BK Sta 1 AC_VI Sta 
2AC_VISta's “ t 


1 AC_VO Sta 

^ 

0 100 200 


± - ^ ^ 

300 400 Time (s') 


Fig. 6 . Injection and/or removal of stations in the WLAN. 
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Fig. 7. Contention window 



Fig. 8. Station throughput 


that the proposed control approach has high accuracy and 
fast convergence speed, and is adaptive to general network 
scenarios. To the best of our knowledge this might be the 
hrst detailed study of using a closed-loop control approach to 
achieve proportional fairness amongst ACs in EDCA WLANs. 
The optimisation of controller parameters is not considered in 
this paper. We leave that for future work. 
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and Af = {1,2, ••• ,n} denotes the set of stations in the 
WLAN. 

As the logarithm of a sum of exponentials is a convex 
function, logV is convex in the transformed variable rj, and 
Ui{ri) is thus concave in r/. ■ 



